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We propose and analyze heralded quantum gates between qubits in optical cavities. They employ 
an auxiliary qubit to report if a successful gate occurred. In this manner, the errors, which would 
have corrupted a deterministic gate, are converted into a non-unity probability of success: once 
successful the gate has a much higher fidelity than a similar deterministic gate. Specifically, we 
describe that a heralded , near-deterministic controlled phase gate (CZ-gate) with the conditional 
error arbitrarily close to zero and the success probability that approaches unity as the cooperativity 
of the system, C, becomes large. Furthermore, we describe an extension to near-deterministic N- 
qubit Toffoli gate with a favorable error scaling. These gates can be directly employed in quantum 
repeater networks to facilitate near-ideal entanglement swapping, thus greatly speeding up the 
entanglement distribution. 

PACS numbers: 03.67.Pp, 03.67.Hk, 03.67.Bg, 32.80.Qk, 42.50.Pq 


Exploiting quantum systems for information process¬ 
ing offers many potential advantages over classical in¬ 
formation processing like highly secure quantum net¬ 
works m and powerful quantum computers [IHS]. One 
of the main challenges for the realization of functional 
quantum computers is to perform gates with sufficiently 
high quality so that the remaining errors can be sup¬ 
pressed by error correction codes, which makes the com¬ 
putation fault tolerant [7]. At the same time, applica¬ 
tions to long distance quantum communication can be 
enabled by quantum repeaters, which combine proba¬ 
bilistic entanglement generation over short distances with 
subsequent entanglement connection steps [3]. Eor these 
protocols, the probabilistic nature of the entanglement 
generation is acceptable, but it is essential that high fi¬ 
delity entanglement is achieved conditioned on a herald¬ 
ing measurement. Experimentally, such high fidelity en¬ 
tanglement is often much easier to implement and may 
be realized in situations where it is impossible to per¬ 
form any quantum operations deterministically. Here we 
introduce a similar concept for gate operations and de¬ 
velop the concept of heralded quantum gates with inte¬ 
grated error detection. In the resulting gate, the infi¬ 
delity, which would be present for a deterministic gate 
is converted into a failure probability, which is heralded 
by an auxiliary atom. Once successful, the resulting gate 
can have an arbitrarily small error. Such heralded gates 
could facilitate fault tolerant quantum computation since 
detectable errors may be easier to correct than unde¬ 
tectable errors HHini- Alternatively it can be directly in¬ 
corporated into quantum repeater architectures for long 
distance quantum communication. 

Optical cavities are ideal for conversion between the 
stationary gate qubits and flying qubits (photons), which 
is fundamental for quantum networks [HHls]. Quantum 


gates can, in principle, also be directly implemented in 
optical cavities [14], but the experimental requirements 
for this are very challenging due to spontaneous emission 
and cavity loss. The essential parameter quantifying this 
is the cooperativity of the atom-cavity system, C. It has 
been argued that directly implementing gates in optical 
cavities leads to a poor error scaling 1 — F oc IjVC, 
where F is the fidelity of the gate HSUS]. However, as 
a result of the integrated error detection, the heralded 
gates that we propose exhibit high fidelities when suc¬ 
cessful. This enables efficient entanglement swapping and 
removes the necessity of intermediate entanglement pu¬ 
rification in quantum repeaters thus increasing the distri¬ 
bution rate significantly. Compared to using other deter¬ 
ministic, cavity based gates, an increase in the rate of up 
to two orders of magnitude can be achieved for modest 
cooperativities (< 100) and a distance of 1000 km [T7] . 


The basic idea is to use a heralding auxiliary atom in 
addition to qubit atoms in the same cavity. One of the 
atomic qubit states, e.g., state |1) couples to the cav¬ 
ity mode while |0) is completely uncoupled (see Eig. la). 
Such a system has previously been considered for two- 
qubit gates Uni [TSHST] . multi-qubit gates [181 122] and 
photon routing [23] . If any of the qubit atoms is in state 
|1) the cavity resonance is shifted compared to the bare 
cavity mode, which can be exploited to make a gate be¬ 
tween two or more qubits by reflecting single photons off 
the cavity m- The efficiency of such schemes, however, 
is limited by photon losses, inefficient detectors and non¬ 
ideal single photon sources [111123]. We circumvent these 
problems by introducing an auxiliary atom in the cav¬ 
ity to serve as both an intra-cavity photon source and a 
detector. As opposed to previous heralded gates in opti¬ 
cal cavities, which relied on the null detection of photons 
leaving the cavity [241126] . the final heralding measure- 
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ment on the atom can then be performed very efficiently. 

In our approach, the auxiliary atom has two metastable 
states l^f), I/), which can be coupled through an excited 
state 1^) (see Fig. pm|). We assume the \E) ^ \ f) tran¬ 
sition to be energetically close to the cavity frequency 
and to be a nearly closed transition, so that we need to 
drive the \g) \E) transition, e.g with a two-photon pro¬ 

cess, (see below). The gate can be understood through 
the phase evolution imposed on the atoms. We con¬ 
sider adiabatic excitation of the auxiliary control atom 
via Stimulated Raman Adiabatic Passage EH EH], driven 
by an external driving pulse with Rabi frequency 
and a coupling to the cavity photon gf. In the case 
when all the qubit atoms are in the non-coupled states 
|00..0), an adiabatic excitation will result in a dark state 
^ —r^|l,/) with zero energy and vanishing phase. 

Here the number refers to the number of cavity photons. 
However, the qubit states T with at least one of the qubit 
atoms in the coupled state, results in a cavity-induced 
shift of the state |1,/, T), which in turn, causes an AC 
Stark shift and dynamical phase to be imprinted into 
the 1^, T) state after the driving pulse is turned off. All 
states but the completely uncoupled qubit state |00...0) 
will thus acquire a phase, the magnitude of which de¬ 
pends on the length of the driving pulse. With an appro¬ 
priate pulse length and simple single qubit rotations, we 
can use this to realize a general Wqubit Toffoli gate or 
a control-phase (CZ) gate. 

Naively, the gates will be limited by errors originating 
from cavity decay and spontaneous emission from the 
atoms, which carry away information about the qubit 
state. These errors are, however, detectable since the 
auxiliary atom will be trapped in state |/) if either a 
cavity excitation or an atomic excitation is lost. Condi¬ 
tioning on detecting the auxiliary atom in state \g) at the 
end of the gate thus rules out the possibility of any dissi¬ 
pative quantum jumps having occurred during the gate. 
As a result, the conditional fidelity of the gate is greatly 
enhanced at the modest cost of a finite but potentially 
low failure probability. 

We now analyze the performance of the gates and de¬ 
rive the success probabilities, gate times and gate errors 
(see Tab. SI in m) The Hamiltonian in a proper rotat¬ 
ing frame is (see Fig. 

H = Ae\E){E\ + gf{d\E){f\ + H.c) + F + 

+ E + 9{d\e)km + H.c), (1) 

k 

where k labels the qubit atoms {h = 1), 2V = Cl\E){g\ 
and we have assumed that all couplings (g,^) are real. 
We have defined Ae = oje — ojg — and Ae = 
uje—ojg—ujL-\- 00 f — ooi^ where cjl is the laser frequency 
and otherwise ujx is the frequency associated with level 
X. We describe the cavity decay and atomic spontaneous 
emission with Lindblad operators so that Lq = y^d cor¬ 
responds to the cavity decay, Lf = fo the 




FIG. 1. (Color online) (a) Level structure of the qubit atoms. 
Only state |1) couples to the cavity and we assume that the 
excited level decays to some level |o), possible identical to 
I/) or | 0 ). (b) Level structure of the auxiliary atom and the 
transitions driven by the weak laser (Q) and the cavity {gf). 
We assume that \E) |/) is a closed transition, i.e. 7 ^ = 0. 


decay of the excited state of the auxiliary atom and 
Lk = y/l\o)k{^\ describes the decay of the excited qubit 
states to some arbitrary ground state |d). The nature of 
|d) is not important for the dynamics of the gates and it 
may or may not coincide with |0) or |1). 

We assume a weak driving pulse justifying for a pertur¬ 
bative treatment of V using the formalism of Ref. [30]. 
In the perturbative description we adiabatically elimi¬ 
nate the coupled excited states of the atoms and the cav¬ 
ity (assuming /Ae Ae and ^ ^), which leads 
to an energy shift of the ground states but otherwise 
conserves them since the Hamiltonian cannot connect 
different unexcited states without decay. The dynam¬ 
ics are therefore described by an effective Hamiltonian, 
HeS = \g){9\ J2n ^nPn where 

A rx [ -g((t-^/2)^+2nC) ] 

and Pn projects on the states with n qubits in state |1). 
For simplicity, we have assumed that the auxiliary atom 
is identical to the qubit atoms such that gf = g and 
Jf = "f (see [29] for a more general treatment) and we 
have defined the cooperativity C = g^ jjK. We consider 
the limit (7^1 and from Eq. § we find that the en¬ 
ergy shift, in the case when all qubit atoms are in |0), 
becomes very small Aq ^ A£;f2^/(167^(7^) ^ 0, i.e., we 
drive into a zero energy dark state as mentioned in the 
description above. On the contrary, for n > 0, the C 
in the nominator of A^ refiects that the coupling of the 
qubit atoms shifts the cavity resonance and as a result an 
AC stark shift of ^ /Ae is introduced. Furthermore, 

we find that in the effective evolution, errors caused by 
spontaneous emission or cavity decay (Lq, Tj, Lk) project 
the system out of the effective space into orthogonal sub¬ 
spaces, which allows for an efficient error detection by 
measuring the ancilla atom. 
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The dynamics described by can be used to im¬ 
plement a Toffoli gate. Assuming the qubit atoms to be 
on resonance (Ae = 0 ) and having Ae ^ 7 a/C gives en¬ 
ergy shifts A^>o ^ r^^/(47V^) while Aq ^ 

Hence, |00...0) is the only state, which remains unshifted 
and we can choose a gate time of tx A'Ky/C^/Vl? to 
make a Toffoli gate. By conditioning on measuring the 
auxiliary atom in state \g) at the end of the gate, the de¬ 
tectable errors from cavity decay and spontaneous emis¬ 
sion only reduce the success probability instead of reduc¬ 
ing the fidelity. Consequently, the fidelity becomes lim¬ 
ited by more subtle, undetectable errors (see Ref. [29]). 
The dominant error originates from the qubit dependent 
decay rate, T^, of ^ |/). As we demonstrate in 
Ref. [29| , this leads to a fidelity lower bounded by 1 — F < 
0.3/C, with a success probability of Pg ^ 1 — 3/v^- 
Thus is a substantial improvement over the leading er¬ 
ror in the case of deterministic cavity-assisted gates. For 
generic states, the fidelity can even be markedly higher, 
and improving with increasing particle number N [29] 

In the special case of only two qubits, the Toffoli gate 
is referred to as a CZ-gate, and in this case, we can even 
improve the gate to have an arbitrarily small error by 
combining it with single qubit rotations. For the gen¬ 
eral Toffoli gate discussed above, we needed Ae = 0 to 
ensure the correct phase evolution, but making the sin¬ 
gle qubit transformations | 0 ) ^ and | 1 ) ^ 

g- 2 (Ai-Ao)t/ 2 |Q^^ at the end of a driving pulse of length 
tcz = k/(A 2 — 2 Ai + Ao)|, ensures the right phase evo¬ 
lution of the CZ-gate without any constraints on Ag. 
Hence, it is possible to tune Ae to eliminate the detri¬ 
mental effect of having a qubit dependent decay rate. 
Choosing Ae = ^\/4C + 1 and Ae = jA e ensures 
Fq = Fi = r 2 , and thus removes all dissipative errors 
from the heralded gate. The conditional error is then 
limited only by non-adiabatic effects, that can in prin¬ 
ciple be made arbitrarily small by reducing the driving 
strength. The success probability is 1 — Pg ^ 6 /v^ in 
the limit C ^ 1 (see Fig. [ 2 a| ). We thus have a heralded 
two qubit gate with arbitrarily small error with a success 
probability that can approach 1 (it is possible to decrease 
the scaling factor of the probability from ^ 6 to ^ 3.4 at 
the expense of an error scaling as l/C by tuning A^;, Ae). 

We now consider the gate time. The gate time of the 
Toffoli gate is tx ^ and for the CZ-gate we 

have tcz ^ 157r\/C7/(2f2^) for C ^ 1. Since tcz > we 
focus on tcz- The gate time is set by the strength (f2) 
of the driving pulse, which is limited by non-adiabatic 
errors. This is investigated in the supplemental mate¬ 
rial where we also verify our analytical results numeri¬ 
cally [29]. Assuming realisitc parameters of /^ = IOO 7 
[23] [3l], we find that a driving of Q = VCjjA keeps 
the non-adiabatic error of the gate below 4 • 10“^ for 
C < 1000. The gate times decreases as Ij^/C as shown 
in Fig. For a cooperativity of 100 the gate time is 
~ 1 /is for typical atomic decay rates. 



FIG. 2 . (Color online) (a) Failure probability (1—Ps - left axis) 
and gate time (tcz -right axis) as a function of the cooperativ¬ 
ity (C) for the CZ gate. The gate time is in units of the inverse 
linewidth I /7 of the qubit atoms. We have assumed a driving 
of H = a/C 7 / 4 . (b) Gate error as a function of the detun¬ 
ing Ae 2 in the two-photon-driven CZ-gate for C = 10, 20, 50, 
and 100. We have assumed that Hmw = 47 ( 7 ^and that 
7 ^ = 7 . The gate error decreases as 7 ^/A |;2 and is indepen¬ 
dent of C. We have assumed H ~ Ae 2 /S resulting in a gate 
time ~ 400 / 7 . Solid/dashed lines are analytical results and 
symbols are numerical simulations (see [29]). For both plots, 
we have assumed k = IOO 7 . 


So far, we have assumed a model where there is no 
decay from |P) ^ \g). In real atoms, there will, however, 
always be some decay |P) -> \g) with a decay rate 7 ^ > 0. 
The result of such an undetectable decay is that both 
the CZ-gate and the Toffoli gate will have an error ^ 
Jg/i'yVC)- To make this error small, it is thus essential 
to suppress the branching ratio 7 ^/ 7 . Below we show how 
to suppress 7 ^ by driving the \g) \E) transition with 

a two photon process. As a result, we realize a CZ gate 
with an error arbitrary close to zero and a Toffoli gate 
with an error scaling as l/C even for a realistic atomic 
system. 

Specifically we think of a level structure for the aux¬ 
iliary atom, shown in Fig. [^ where we still assume 
|P) ^ I/) to be a closed transition. For simplicity, we 
have also assumed IP 2 ) ^ \g) to be a closed transition. 
Such a level structure could, e.g. be realized in ^^Rb as 
shown in Fig. [^ We assume that a microwave field cou¬ 
ples the two excited states such that we can have a two 
photon transition from \g) \E) and that O is small, al¬ 

lowing for a perturbative treatment of the coupling. Thus 
we can map the system to a simple three-level atom with 
levels 1^), |P) and |/) and a decay rate 7 ^ and drive 
between \g) and |P), determined by the two photon driv¬ 
ing process as shown in Fig. The dynamics are thus 
similar to what we have already described for the sim¬ 
ple three level atom except that we have the extra decay 
7 ^ that introduces an error in the gates ^ (%/ 7 )/v^, 
as previously described. In the limit C ^ 1, we find 
Igh ^ 4 !/^^ • Thus by increasing Ae 2 ^ we can in prin¬ 
ciple make these errors arbitrarily small. The error of 
the CZ-gate for different Ae 2 is shown in Fig. |2b[ as- 
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FIG. 3. (Color online) (a) Level structure of the auxiliary 
atom and the transitions driven by a weak laser (Q), a mi¬ 
crowave held (Qmw) and the cavity (g/). We assume that 
\E) o I/) is a closed transition and, for simplicity, we also as¬ 
sume that \E 2 ) o I^) is a closed transition but this is not a ne¬ 
cessity. Here with r = 1, 2, 3 refers to how the atom 

may be realized in the (5^P3/2) states \E^^^ = = r) 

5 ^5 'i/2 of Rb®^. (b) Effective three-level atom realized by 

mapping the two-photon drive to give an effective decay rate 
7 ^ and an effective drive Q. 


suming an initial state of (|0) + Note that in 

order to prevent an increasing scattering probability of 
level |P 2 ), we need to have ^mw oc (7^/^ resulting in a 
gate error that is independent of the cooperativity. m- 
The success probability and time of the gates are the 
same as before with Q ^ Q ^ With similar 

considerations about the validity of our perturbation as 
before, we hnd that for realistic parameters, we can use 
Q = A£; 2 / 8 ,Omw ^ resulting in a gate time of 

^ 10 /rs for typical atomic decay rates and C < 1000 [29] . 

As an example implementation, we consider ultra-cold 
^^Rb atoms coupled to nanophotonic cavities [23[ [3T] . 
There are some additional errors originating from the 
extra states in the ^^Rb atoms in this case. In Ref. [29] . 
we treat these errors and find that with a detuning of 
^E 2 = IOO 7 and a cooperativity of (7 ~ 100, a heralded 
CZ gate with ^ 67% success probability and a heralded 
error of ~ 10“^ can be realized in 10 /is time. This 
justifies neglecting atomic decoherence which is typically 
much slower. Alternatively the gate can be implemented 
with atom-like solid-state qubits such as NV and SiV cen¬ 
ters in diamond [32]. These systems can exhibit closed 
transitions and long-lived electronic spin states which are 
the essential requirement for the gate [33] , while high co- 
operativities are possible in diamond nanocavities [34] . 
A particular advantage of such system is the long-lived 
nuclear spin degrees of freedom, which allows each of the 
color centers to act as a multi-qubit quantum network 
node [ 35 ]. By entangling electronic spins via the her¬ 
alded gate, a high-fidelity, fully deterministic gate can 
subsequently be performed on qubits stored in nuclear 
spins m- 

As a particular application, we consider a quantum 
repeater where entanglement is first created in small seg¬ 


ments (links), which are subsequently connected using 
entanglement swapping [36]. By organizing the repeater 
in a tree structure, the probabilistic nature of the gate 
can be efficiently circumvented. The success rate of dis¬ 
tributing entanglement across the total distance L, scales 
as ^ (I//I/o)^“^°^ 2 ( 3 /p)^ where p < 1 is the success prob¬ 
ability of the swap, L is the total distribution distance 
and I/O is the length between the links [17] (note that in 
the limit p ^ 1 , the above expression underestimates the 
rate, e.g., for p = 1 the actual rate is ^ 3 times faster for 
128 links). This is a substantial improvement over direct 
transmission where the success rate scales exponentially 
with L. For a realization with nuclear spin memories 
where the swap can be performed deterministically the 
rate can scale even better as ^ \og 2 {L/Lq)~^. In order to 
maintain the favorable scaling without resorting to time 
consuming purification, the total number of links, Amax 
should be kept below A%ax -ln(Ffinai)/(eo+e 3 ), where 
Ifinai is the required fidelity of the final distributed pair 
and eo, ^ 1 are the errors of the initial entanglement 
generation and the entanglement swapping respectively. 
Thus, it is essential that the errors are kept small, which 
can be obtained with the heralded gate. 

In conclusion, we have introduced a heralded two-qubit 
quantum gate with a conditional fidelity arbitrarily close 
to unity and an A^-qubit Toffoli gate with an error scaling 
as 1/(7. The gates have a built-in error detection process, 
which removes the necessity of extracting the error by the 
more complicated process of entanglement purification or 
quantum error correction. Our gate is designed for the 
specific case of optical cavities, and allows exploiting re¬ 
alistic systems for quantum communication, even though 
the error rate would inhibit this with deterministic gates. 
Similar advantages can be realised in other systems such 
as those based on circuit QED, where certain errors could 
be heralded and thus alleviate the daunting requirements 
of fault tolerant computation. 
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SUPPLEMENTAL MATERIAL 


This supplemental material to the article ’’Heralded quantum gates with integrated error detection in optical 
cavities” describes the details of the perturbation theory and the derivation of the effective Hamiltonian ^eff and 
effective Lindblad operators. We describe the situation both with and without a two-photon drive. Furthermore, 
we present the results of a numerical simulation of the full dynamics of the gates to verify the results found with 
perturbation theory and address the question of how strong a drive we can allow for. In the end, this determines 
the gate time as described in the article. Finally we discuss of the additional errors described in the final part of the 
article. 


PERTURBATION THEORY 


We will now give the details of the perturbation theory and the derivation of the effective 
the success probabilitites, gate times and gate errors (see Tab. SI). Our perturbation theory 
operator formalism described in Ref. m- 


operators together with 
is based on the effective 


Gate 

Origin of error 

Error 

Probability 

Time 

GZ-gate 

79 = 0 

79 > 0 

0 

'-yVc 


Ibn'/C'y 
^ 202 

Toffoli 

r* ^ r,- 
79 > 0 

< 0.3 
~ C 

^Vc 

- 1 - ^ 

47r7C7 
^ 02 


TABLE SI. The errors, success probabilities and gate times of the iV-qubit Toffoli gate and the CZ-gaie considered in the 
article. Note that the branching fraction 7^/7 can be made arbitrarily small using a far detuned two-photon driving as explained 
in below. T^ is the rate of detectable errors for the qubit state with i qubits in state |1). The success probability of the CZ-gate 
can be increased at the expense of an error scaling of l/C as explained in the article. 


First, we treat the simplest situation where the auxiliary atom is directly driven to an excited state \E) by a weak 
classical drive U as shown in Fig. (reproduced from Fig. 1 in the article). Note that we allow for some decay from 
\E) \g) with decay rate 7 ^ as opposed to the situation in the article. We will later consider the situation where this 

decay rate is suppressed using a two-photon drive. The level structure of the qubit atoms are also shown in Fig. [S^ 



FIG. S4. (a) Level structure of the qubit atoms. Only state |1) couples to the cavity and we assume that the excited level 
decays to some level |o), possible identical to |/) or |0). (b) Level structure of the auxiliary atom and the transitions driven by 
the weak laser (Q) and the cavity (gj). We allow for some decay from \E) \g) with decay rate 7 ^. 


The Hamilton describing the system in a proper rotating frame is given by Eqs. (l)-(3) in the article and is 
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reproduced here 


H = He + V + V\ 

He = AE\E){E\+gf{a\E){f\ + H.c) 

+ Ae\e)k{E\ + 5 (d|e)j;(l| +H.c), 

k 

V=^\E){9\, 


(53) 

(54) 

(55) 


where we have assumed for simplicty that all couplings (^, are real and k labels the qubit atoms {h = 1). We have 
defined Ae = uje — and Ae = uJe—^g—— where ujl is the laser frequency and otherwise is 

the frequency associated with level x. Note that we assume the cavity frequency to he ujc = ool ^g — such that 
we are on resonance with the \g) \E) \ f) two-photon transition. 

The dissipation in the system is assumed to be described by Lindblad operators such that Lq = y/kd describes the 
cavity decay with decay rate Lg = y/^\g){E\, and Lf = y^\f){E\ describes the decay of the auxiliary atom and 
= ^/l\o)i{e\ describes the decay of the qubit atoms {k = 1,2.. .N). As described in the article |d) may or may 
not coincide with | 0 ) or | 1 ). Assuming that ft is weak (p?/A e Ae and ft g), we can treat the driving as a 
perturbation to the system. As shown in Ref. m, the dynamics of the system is then governed by an effective master 
equation of the form 


p = i 


P, ^eff 


+ Lfp{Lf)^ - 1 ({Lf)^Lfp + p{Lf)^Lf) , 


(S 6 ) 


where p is the density matrix of the system, ^eff is an effective Hamiltonian, and are effective Lindblad operators 
with X = 0, ^, /, /c. The effective operators are found from: 


ffeff = 


1 

NH 


(-^nh)^) 


V 


Lf = 


where 




(57) 

(5 8 ) 

(59) 


is the no-jump Hamiltonian. The Hilbert space of the effective operators can be described in the basis of {|^), |/)} 
of the auxiliary atom and the states {|0), |1), |d)} of the qubit atoms. To ease the notation, we define the projection 
operators Pn which projects on to the states with n qubits in state | 1 ). From Eq. (S7) and (S 8 ) we then find: 


N 


iAe /2 + nC 


^-ft\ 
iLeff = / -Re \ 

^=0 \Ae{iAE/2 P Cf) E AetiC ^ 


\g){g\^Pn 


N 


^ ^ 1 ^) (^1 ^ Pn 


(SIO) 


n=0 

N 




^/CfAeQ 


n'^0 Ae{iAE/2 + Cf) + uAeC 

N 

E<n\f){9\®Pn 

n=0 
N , 


\f){9\®Pn 


(Sll) 


^ eff 


{iAe/2EnC)ft yfTg 
ApAEp E Cf') E uAeC 7 

N 


= E 


\g){g\®Pn 


= EEn\g)ig\®Pn 


(S 12 ) 


^ eff 


n=0 

^ 1 

= E 


(zAe /2 + nC)n fyj 
2 Ae(zA£:/2 + C/) + nA^jC 7 


\f){g\®Pn 











(S13) 


N 




n=0 

N 


Lf = y-^ 

^ 9 . / 




^ 2.^/7 Ae(iA£:/2 + C/) + nA£:C 

N 

= Y.^n\f){9\®\0)k{l\®Pn, 

n=l 




(S14) 


where we have defined the cooperativities (7(j) = for the qubit (auxiliary) atoms and the complex detunings 

A £;7 = A^; — ijfl2 and Ae 7 = Ag — '^ 7 / 2 . Note that we have defined the parameters ^o!n 5 Elqs. 

) to characterize the decays described by the Lindblad operators. Note that = 0. In our calculations 
we parametrize the difference between the auxiliary atom and the qubit atoms by C/ = aC and 7 / = to easier 
treat the limit of C ^ 1 that we are interested in. 


(S11)-(S14 


Success probability and fidelity 


Eqs. (Sll )-( [ST4| ) show that the effect of all Lindblad operators, except is that the state of the auxiliary atom 
is left in state |/). All these errors are thus detectable by measuring the state of the auxiliary atom at the end of the 
gate. For the heralded gates where we condition on measuring the auxiliary atom in state \g) at the end of the gates, 
these detectable decays therefore do not effect the fidelity of the gates but only the success probability. The rate T^ 

and assuming an initial 


+ 


of the detectable decays for a state with n qubits in state | 1 ) is T^ 
qubit state described by density matrix pqubit the success probability of the gates is 


' 0,n I 


^eff 

f,n 


+ r: 


,eff I 


-^success — ^ Tr {e-r"*^-Vqubit7} , (S15) 

n=0 

where tgate is the gate time and Tr denotes the trace. 

Having removed the detrimental effect of the detectable errors by heralding on a measurement of the auxiliary atom 
the fidelity of the gates will be determined by more subtle, undetectable errors (see below). We define the fidelity, F 
of the gate as 

F=j^mg\p^uut\gm, (S16) 

-^success 

where we have assumed that the ideal qubit state after the gate is a pure state and pqubit is the actual density 
matrix of the qubits and the auxiliary atom after the gate operation. 


A-qubit TofFoli gate 


As shown in the article, the effective Hamiltonian in Eq. (SIO) is sufficient to make a Toffoli gate by putting the 
qubit atoms on resonance (Ag = 0). We will now treat the worst case and average fidelities of the general Toffoli gate 
referred to in the article. The undetectable errors limiting the fidelities are the following. 


• As described in the article the energy shifts of the coupled qubit states are all A^>o ~ /{4J^/C) in the 

limit C ^ 1. However, to higher order in C, we find corrections on the order to the energy shifts, 

which depend on the number of qubits that couples. The gate time of the Toffoli gate is tx r\j FK\fC ^and 
consequently, the higher order corrections give uneven phase shifts on the order of 0{C~^) for the coupled qubit 
states at the end of the gate. This leads to a phase error in the fidelity of 0{C~‘^). 

• The difference between the rates of detectable errors (T^) for different qubit states changes the relative weight 
of the qubit states during the gate. This error wil be 0{C~^) as shown below. 

• For 7 ^ > 0 the undetectable decay from \E) \g) in the auxiliary atom will destroy the coherence between the 

qubit states. We find that this error will be ^ now, we will assume that 7 ^ = 0 and thus ignore this 
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FIG. S5. (Color online) (a) Gate error of the Toffoli for different initial states as a function of cooperativity. We have plotted 
the generic error for N = 5,10, and 15 and the upper bound of the error. Note that the generic error decreases as N increases. 
We have hxed A such that Fo = Fi and have assumed that a = = 1. (b) The failure probabilities 1 — Psuccess,up and 

1 — iTuccess,gen as a functiou of cooperativity. 1 — Psuccess,gen is plotted for iV = 5,10,15. We have used the same assumptions 
as in (a). In general, the failure probability only have a weak dependence on N. Note that the line for 1 — Psuccess, gen,N = 5 

coincides with 1 Psuccess,up* 


error since we will show that we can suppress the branching fraction 7^/7 arbitrary close to zero by having a 
two photon driving. 


Assuming that 7 ^ = 0, the dominating source of error limiting the performance of the Toffoli gate is thus the 
difference between the rates of the detectable errors for the qubit states. We tune Ae such that Fq = Fi and the 
largest difference between the detectable errors is thus between the completely uncoupled state and the state with all 
qubit atoms in state |1). As a result, we can find an upper bound on the fidelity of the N qubit Toffoli gate, considering 
an initial state |0)^^ + the limit N ^ 00 because this state experiences the largest difference between the 

number of coupled and uncoupled qubits. We find that the upper bound on the fidelity and the corresponding success 
probability is 


16(q^ T / 3 ) C 
{a + 2i3)7t 1 

2\/SV O' /3 \/C 


(517) 

(518) 


. In general, the fidelity of the gate will, however, be larger than what is suggested above. Considering a generic input 
state (| 0 ) + with the same parameters as above, we find 


1 - k{N) 


1 


■ success,gen 


1 - 


{d{N)a + 2 / 3 ) 7 r 1 
2v^\/ O' P ^/C 


(519) 

(5 20 ) 


where k{N),d{N) are scaling factors which depend on the number of qubits N. We calculate k{N) and d{N) 
numerically for A" = 1 — 100 using the perturbation theory and find that that they both decrease with N (see Fig.|S5[). 
The upper bounded and generic fidelities and corresponding success probabilities are shown in Fig. for different 
number of qubits, N. As N increases we obtain higher generic fidelity, whereas the success probability is almost 
independent of N. 


CZ-gate 

In the special case of only two qubits the Toffoli gate is referred to as a control-phase (CZ) gate. As shown in 
the article, we can, in this case, completely remove the errors from the gate by choosing the detunings Ae and Ag 
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such that ro = ri=r 2 and combining it with single qubit rotations we can ensure the right phase evolution. In the 
general case where 7 ^ 1 , the detunings Ae and Ae are 


Ae 

Ae 


aCj‘^ 

2Ae 


(5 21 ) 

(5 22 ) 


The success probability of the gate is then 


1 — TT 


^/c' 


(S23) 


and we find that the gate time is tcz — VC in the limit C ^ 1 . 


Two-photon driving 


We now describe the details of the implementation where the auxiliary atoms is driven by a two-photon process 
as shown in Fig. S 6 (reproduced from Fig. 4(a) in the article) in order to suppress the dominant undetectable error 
caused by spontaneous decay of the auxilliary atom into the state \g) {Lg). 




(a) |g) = |1,1) (b) |g) 


FIG. S 6 . (a) Level structure of the auxiliary atom and the transitions driven by a weak laser (Q), a microwave field (Qmw) 
and the cavity (gj). We assume that \E) |/) is a closed transition and for simplicity we also assume that IF 2 ) l^') is a 

closed transition but this is not a necessity. The figure also indicates how the levels could be realized in ®^Rb. Here 
with r = 1,2,3 refers to state |F^^^ = r, mp^^ = r) in 5^*S'i/2 (5^P3/2). (b) Effective three level atom realized by mapping the 
two-photon drive to an effective decay rate 7 ^ and an effective drive Q 

The Hamiltonian in a proper rotating frame is 

H = He + V + v\ (S24) 

He = 1 ^e\E){E\ + 11e2\E2){E2\ + 9f{a\E){f\ + H.c) 

+ ^{\E){E2\+ H.C.) 

+ ^Ae|e)fe(e| + 5 (a|e)fc(l| +H.c.), (S25) 

k 

V=^\E2){9\, (S26) 

where we have now defined Ae = LOE - - l^laser - WmW, Ae 2 = i^E2 - Wg - OJlaser and Ae = UJe - Ojg - UJlaser - 

bJjviw + bJ/ — wi. Here coiaser IS the frequency of the laser drive (f2), wmw is the frequency of the microwave field 
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(^Mw) and otherwise ujx is the frequency associated with level x. We assume that the frequency of the cavity is 
= ^laser + <^MW ^ ~ ^uch that the thrcc photon Raman transition from \g) \f) is resonant. We have 

assumed that Ae 2 is large and positive such that the rotating wave approximation is valid for the microwave field. 
The Lindblad operators describing the system are the same as described below Eq. ( |S3| ) except that Lg y^\g){E 2 \. 
Assuming a weak drive 11, we can follow the same recipe as before to find the following effective operators describing 
the dynamics of the system. 


N 






Ae{iAE/‘2 + Cf) + uAeC 




^ 47 [A£:2Ae(zA£:/2 + C/) + nA£:AB 2 C-zAell^w/ 8 -«f^MW ^/4 

= E^n^\9){9\®Pn 

— 1 ^/Cf 


\g){g\®Pn 


n=0 

N 


i-r '=E 


^ AE2Ae{iAE/2 + Cf) + nAE2AEC — iAgO^w/^ “ 

N 

= E'^tn\f){9\®Pn 

H Ae(iA£;/2 + + 7iA£;(7 


\f){9\®Pn 


n=0 

N 




_Q 2 A£:2Ag (iA£:/2 + )-j-?2A£:2A^jC^ iAg^^^/8 7 


\g){9\<2)Pn 


N 


= J2^f!n^^\3){9\<^Pn 


n=0 

^ O 

^eff(2) _ 

/ A-/ 4 


(iAe /2 + n(7)llMW 




^ 4 AE2Ae{iAE/2 + Cf) + nAE2AEC — iAell^w/^ “ 


|/)(^| ( 8 ) Ai 


AT 


= E"^fn\f)i9\®Pn 


n=0 

AT-l 


H M W ^ 


^eff(2) ^ _ 

^ 4^ AE2Ae(zAE/2 + C/) + nAE2^EC - iAettl^y^/8 - nQl^y^C/4: 

N-l 


(S27) 


(S28) 


(S29) 


(S30) 


\f){9\ ® |o)fe(l| ® P„ 




(S31) 


n=0 


eff(2) 


' 9,ri 


eff(2) j eff(2) 

and rn 


where we have defined the complex detuning Ae 2 J = Ae 2 — and the parameters 

to characterize the decay described by the Lindblad operators. 

We are interested in the limit of large detuning Ae 2 and large cooperativity C. In this limit, we find that the 
dynamics of the system can be mapped to a simple three level atom with effective driving Q ^ QQmw/{‘^^E 2 ) and an 
effective decay 7 ^ 


^ 7 ^^mw/^I ;2 shown in Fig. S 6 In principle, the effective operator leads to an effective 

decay rate of 7 = 7 ^ 11 ^/A |;2 to lowest order in C but we find that this first order term do not destroy the coherence 
between the qubit states since it is independent of n. There is, therefore, no effect of these scattering events and the 
performance of the gate behaves as if there is an effective decay rate of 7 ^ ^ 7 ^^mw/^I; 2 - Note that we also have an 
AC stark shift imposed on the level \g) by the laser characterized by H. This will give an overall phase to the system, 
^ ( 4 A£; 2 )t, which we can neglect since it does not influence the gates. Since we can do the mapping to the simple 

three level atom, we find similar results for the performance of the gates for the two-photon scheme as for the simple 
three level scheme only with effective decay 7 ^ and drive given by the two-photon process. Note, however, that we 
now assume 7 ^ > 0, which introduces an undetectable error as previously mentioned. We find that this introduces an 
error in the fidelity of both gates of roughly 

{a‘^ — AaP — 6yd^)7r^ 7 ^ 


128/32 


__ + 4(4/3 + 6P^)7r Jg^MW 1 

l‘^^E2 16v^(a + 2 /?)(Q; + 5/?)) 7 AI 2 VC' 


(S32) 


Nonetheless, this error can be suppressed arbitrarily much be increasing Ae 2 ^ which enable us to have a heralded 
CZ-gate with arbitrarily small error in a realistic atomic setup using the two-photon drive. 
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GATE TIME 


Here we address the question of how strongly we can drive the system and still maintain the validity of perturbation 
theory. We need to adress this question since the gate time depends inversely on the driving strength as shown in the 
article and hence this limits the achievable gate time. A necessary criterion for our pertubation theory to be valid is 
that the energy shifts (see Eq. ( |S1Q[ ) and Eq. ( |S27[ )) are small compared to the driving, i.e. ^ ^ 1 . Erom 

Eq. (SIO) we find that A^/0^ ^ 12^/(16A|.) to leading order in the cooperativity C and this criterion is therefore 
met for <C dA^;. Similarly, for the two-photon process, we find from Eq. (S27) that r\j 11 ^/ ( 16 A|. 2 ), to 

leading order in C. Here we thus need 12 ^ 


Another criterion need to be meet in order for our perturbative theory to be valid. If none of the qubits couple, we 
are effectively driving the auxiliary atom/cavity system into a dark state of the form cos(^)| 0 ,^) — sin(^)|l, /), where 
the mixing angle is 6 ^ Q./g. Here the number refers to the number of cavity photons. To adiabatically eliminate the 
state I/, 1) from the Hamiltonian as we have done in the effective Hamiltonian requires that Vt/g 1 . Mapping this 
criterion to the effective three level scheme realized in the two-photon scheme gives OOmw/( 2 A£; 2 ^) ^ 1 . 

Einally, we need to consider the scattering of photons from the level \E2) in the two-photon scheme. If the number 
of scattering events, riscat is large compared to 0 ^/(I 6 A|. 2 ), the perturbation theory is not valid even though the 

other criterions are met. We find that riscat ^ ^ for fho CZ-gate and we thus need to have A 02 — ^ 1 

The different criterions for the validity of the perturbation theory are summarized in Tab. |S 2 | 


Simple scheme 

Two-photon scheme 

0/(4Ab) < 1 

12 /(4A£;2) ^ 1 

Q/g 1 

1212 mw/(A£;25') 1 

- 

3'\/U7^12^/(A|;212mw) ^ f 


TABLE S2. The criterions for our perturbation theory to be valid. 


Eor all the gate schemes, we have assumed that Ae oc VC^. The first criterion for the simple scheme (see Tab. S 2 ) 
can thus be met with a driving of 12 = ajVC where a/4 ^ I. This driving results in a gate time that decreases 
as The value of a will determine the size of the non-adiabatic error. The second criterion is also fulfilled for 

this driving as long as ^ I. In realistic systems such as the nanocavity system described in Ref. [23l [31], the 

ratio can be on the order of 100-1000. Eor the two-photon scheme, we find from Tab. S 2 that we can choose 
12 = a2A£;2,12MW = 67 ( 7 ^/^ where the constants a 2,6 will determine the size of the non-adiabatic errors as before. 
Similar to the situation in the simple scheme we assume that 1 . 


NUMERICAL SIMULATION 


In order to confirm our results, we numerically integrated the full Master equation, defined by the Hamiltonian H 
and the Lindblad operators, Lj G {Lq, Ej, Li, I/ 2 }, 


j^p{t) = ^ i [2Ljp{t)L] - p{t)L]Lj - L]Ljp{t) 


(S33) 


We used the QuTiP 2 package m, for Python, to set up the problem and used its 12th-order numerical integration 
algorithm to find the solution p{t) as a time series. Then, we used the routines of the same package to analyze the 
results (see Ref. EH). 

Eor each time series p(t), we determined the gate time tgate? fhe success probability Psuccess? and the fidelity F. We 
picked |' 0 o)i 2 = 172 ( 1 ^^ ^ l^))i ^ 171 ^ l ^))2 initial state of the two qubits, \g) for the control atom, and 

zero photons in the cavity. Starting from here, we let the system evolve under the Master equation Eq. (S33), and 
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determined Pg{t)^ the conditional state pg{t) and F{t) as a function of time: 


Pg{t) =Tr p{t)\g){g\ 


Pg(t) = 


\9){9\p{t)\g){g\ 


Pgit) 


Fit) = max 
01,02 




, 01,02 


Tic, c {^Pg{^)) 




01,02 


where Trc, c is the partial trace operation over the control atom and the cavity, and is the target state 

transformed with two single qubit ^^-rotations: 

= c/i(<^i);72((/.2)i(|oo) + |oi) + |io) - 111 )), 

where Uk{4^k) = exp [i|l)/c(l|/c0/c] is the 2 :-rotation of qubit k {= 1,2) by the angle (j)k- From these time series, we 
determined the gate time tgate by finding the timepoint where F{t) is maximal, 

^gate = argmax F{t) 

t 


The fidelity and the success probability of the gate is then defined as F = ^(tgate), F^success = F^^(^gate)- 

Plots of Fig. show the gate time (tgate) and the success probability (Psuccess) as a function of a, for 7 ^ = 0, 
7 = O.Ol/D, O = a^\/C with C G {10,30,100,300,1000}. The detunings, A^; and Ae were chosen to be close to 
their optimal value, determined from the adiabatic theory, and numerically optimized to result in identical effective 
1^) ^ I/) transition rates Fq = ri = r 2 for the qubit sectors | 00 ), | 01 ), | 11 ). The rates Vj were found by numerically 
diagonalizing the master equation for the qubit sectors separately, and finding the eigenvalue with the smallest (but 
non-zero) absolute real part. This numerical optimization yielded the maximal fidelity. The symbols correspond to 
the numerical result, whereas the solid lines show the theoretical values. The agreement of the results confirms the 
validity of the adiabatic theory for a driving Q. = a^^/C for C < 1000 and a < 0.25. Note, however, that Fig. S7 shows 


how the success probability deviates from the adiabatic result for a > 0.25. A weak increase of this deviation with 
the cooperativity is seen but from simulations at high C, we believe that this can be removed by gradually ramping 
Vt up and down to maintain adiabaticity at the beginning and end of the driving pulse. This was not included in the 
simulations behind Fig. for simplicity. 




FIG. S7. Gate time (left) and failure probability (right) as a function of driving strength (a) for — 0.01, 7 ^ = 0, and 
C G {10, 30,100, 300,1000}. The driving strength was assumed to be = a^\fC. 


Fig. shows the conditional mfidelity of the gate as a function of a for the same parameters. The simulation 
confirms that using a = 0.25 is enough to push the (conditional) infidelity of the gate below 4 • 10“^. The fidelity 
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is limited by non-adiabatic effects, which can be suppressed by decreasing Q as shown in the figure. Adiabatically 
ramping up and down at the beginning and end of the gate will also improve the adiabaticity but, for simplicity, we 
have not included this in the simulations described here. Note, however that for high C {C > 1000 ), we find that this 
gradual ramping of Q significantly decreases the non-adiabatic error. 



FIG. S 8 . Conditional infidelity of the gate as a function of the driving strength (a) for — 0.01, 7 ^ = 0, and C G 
{10,30,100,300,1000}. The shaded region (at ~ 10“®) shows the limit of numerical accuracy. The driving strength was 
assumed to be = ajVC. 

We repeated the above analysis for the two-photon-driving Hamiltonian in Eq. ( |S24| ). We chose "f = Jg = 7 / = 
0.01 / 4 , and Hmw = and chose Ae and Ae detunings again close to their adiabatic optimum, but 

numerically optimized them with the same procedure as previously. Plots of Fig|S9|show the gate time and the success 
probability as a function of Ae 2 for C G {10,20,50,100}. Symbols indicate the numerical results while solid lines 
show the theoretical values. Fig. 2b in the article shows the conditional infidelity of the two-photon-driven gate as a 




FIG. S9. Gate time (left) and failure probability (right) as a function of Ae 2 for 7 = 7^ = 7^ = O.Oln, Q = 
Qmw = 47C^/^ C = { 10 , 20 , 50,100}. 


function of C for the same choice of parameters. With these results, we confirm that by increasing Ae 2 we can lower 
the infidelity error to an arbitrary small level. The gate time is a constant of the cooperativity since we have increased 
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^MW as in the simulations. We do see some deviation from the analytical results due to non-adiabatic effects, 
which could be suppressed by decreasing ft at the expense of an increase in the gate time. Finally a gradual ramping 
of ft could also decrease the non-adiabatic errors but for simplicity, we have not included this in the simulations. 


ADDITIONAL ERRORS 

There are some additional errors in a realistic atomic setup that we have not treated in detail so far. Here we 
estimate the dominant errors and determine under which conditions, they can be sufficiently suppressed such that 
they do not limit the performance of the gates. We assume a realistic atomic setup where ^^Rb atoms are used 

both for the auxiliary atom and the qubit atoms. In the ^^Rb atoms, we assume that \g) = |1,1), |/) = |2,2) and 

\E 2 ) = |2®,2®), \E) = |3®,3®) where with r = 1,2,3 refers to state = r, mp^ = r) in 5^5'i/2 (5^P3/2). 

In this case, we estimate that the dominant errors are: 

• In our perturbative theory, we have assumed that the laser field {ft) only couple \g) \E2) in the auxiliary 

atom. However, for a large detuning Ae 2 it may also couple |/) ^ |P), which could lead to an undetectable 

error where the auxiliary atom is pumped back to \E2) from, which it decays to \g). This error is, however, 
suppressed by the large frequency separation, of \g) and |/), which is ^ IOOO 7 for Rb. We estimate 
the error using effective operators to find the decay rate back to \g), assuming that the auxiliary atom starts 
in \E2) and treating the drive as a perturbation while neglecting the cavity coupling. This is valid as long 
as Ag ^ Ae^ which is fulfilled for C < 10000 since Ae ^ VC. The error increases with Ae 2 but even for 
Ae 2 ^ 4 OO 7 we find that for Omw = 47(7^/^, = Ae 2 /^ the error is < 10“^. 

• The microwave might also couple the ground states |0) — |1) of the qubit atoms and the ground states \g) — \f) of 
the auxiliary atom. The coupling of |0) — |1) means that the qubit atoms also couple to the cavity even though 
they are in state |0). We estimate the error from this to be on the order of ^mw/ i^g — {^E 2 — E A 2 ^ 3 ))^ 
where A 2^3 is the splitting between IP 2 ) and \E). For ^^Rb, A 2^3 ~ 447 . Below we argue that we need 
A^; < 0. Since Ae ^ —\fC^ this error will increase slowly with cooperativity but it is suppressed by Ag. For 
^MW = 47 ( 7 ^/^, we find that the error is < 10“^ for C < 1000 even for Ae 2 ^ 4 OO 7 . The errors from the 
coupling of the states \g) — \ f) in the auxiliary atom will likewise be suppressed by the large energy splitting 
A^. These errors can also be further suppressed by decreasing ftuw at the cost of a larger gate time. 

The above errors can be highly suppressed using e.g.^^Sr, ^^^Ba+ or ^^Ca+ instead of ^^Rb. For these atoms, the 
ground states can be encoded in the Sq and Pq manifoldes for ^^Sr and the S 1/2 and 1 ^ 3/5 manifolds for ^^^Ba+ and 
^^Ca+, which have separations at optical frequencies between the stable states. 

A final error that we will consider is that the transition \E) |/) will not be completely closed if the cavity is 

linearly polarized. This will, e.g. be the case for the system in Ref. [31]. Such a cavity also couples |/) to the 
states |1®, 1®), |2®, 1®) and |3®, 1^). From |1®, 1®) and |2®, 1®) there might be an undetectable decay back to \g), which 
will introduce an error oc If^/C in the gates. The probability of an undetectable decay from these states should be 
compared to the probability of the detectable decays where the cavity photon is scattered of the qubit atoms instead. 
For ^^Rb, we estimate this error by compairing the strengths of the effective couplings from |/) to |1®, 1®) and |2®, 1®) 
with a subsequent decay to \g) with the strength of the effective coupling from |1) to \e) in the qubit atoms with a 
subsequent decay back to |1). The latter process has a detuning of A^; while the first two are additional detuned by 
the energy gaps between |3®, 3®) and |1®, 1®) and |3®, 3®) and |2®, 1®) respectively, assuming that A^; < 0. We find that 
since | A^;! grows as \/C the error increases from ^ 5 • 10“^ at (7 = 1 to a maximum value of ~ 2 • 10“^ for C ^ 3000 for 
which A^; is comparable to the extra detunings of the |1®, 1®) and |2®, 1®) transitions compared to the |e) transition. 
For C > 3000 the error decreases as IjVC. Note that this error could be removed by making a 4 photon drive from 
1^) to \E) by letting \g) = |1, — 1). Another approach is to consider other atoms such as ^^Ca+, with more favorable 
levelstructures. The state \g) could be encoded in the 3^D^j2 subspace while the state |/) could be encoded in the 
4^S'i/2 subspace and similarly for the qubit states |0) and |1). In such a setup, we will have separations of optical 
frequencies between the qubit states and we can remove the decay from the excited state back to \g) by, e.g. driving 
from 3^D^j2 to 4 ^Pi /2 through 3^D^j2. 


